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Abstract. We consider the problem of enumerating planar constellations with two points 
at a prescribed distance. Our approach relies on a combinatorial correspondence between 
this family of constellations and the simpler family of rooted constellations, which we may 
formulate algebraically in terms of multicontinued fractions and generalized Hankel determi- 
nants. As an application, we provide a combinatorial derivation of the generating function of 
Eulerian triangulations with two points at a prescribed distance. 
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1. Introduction 



From the initial work of Hurwitz ( 1891 ) about transitive ordered factorizations of the identity 



in the symmetric group, to the bijective approach of Bousquet-Melou and Schaeffer (2000) and 



Bouttier et al. (2004), including the more algebraic approach of Goulden and Jackson (1997), 



constellations appear in many forms in different areas of combinatorics. We refer the reader to 



the book of Lando and Zvonkin ( 2004 ) for an extensive review of the variety of contexts in which 



constellations appear. 

In this paper, we focus on the map point of view, and consider the problem of enumerating 
planar constellations with two points at a prescribed distance. In Bouttier et al. (2004), this 



problem has already been tackled by giving a bijection between this family of constellations and 
a family of decorated trees called mobiles, which results in recurrence equations that characterize 
the associate generating functions. Remarkably, these recurrence equations admit explicit solu- 



tions (Di Francesco 2005 Section 6.7), reminiscent of integrable systems (see also Bouttier et al. 
(2003a) for the solutions of the equations corresponding to 2-constellations). The motivation 



of our work is to explain combinatorial!}' the form of these solutions. Following the approach 



of Bouttier and Guitter (2011), we exhibit a connection between our enumeration problem and 
the a priori simpler problem of counting rooted constellations. This connection relies on a de- 
composition involving some lattice paths, which we may rephrase in terms of "multicontinued 
fractions", in the spirit of Flajolet (1980). Solving our enumeration problem amounts to finding 
unknown coefficients in such a multicontinued fraction, which can be done via a generalization of 
Hankel determinants. We illustrate this approach in the case of Eulerian triangulations, where 



we provide a combinatorial derivation of a formula found in Bouttier et al. (2003b) for the gen- 



erating function of Eulerian triangulations with two points at prescribed distance. These objects 



are in bijection with "very- well-labeled trees with small labels", see also (Bousquet-Melou, 2006 



Section 2.1). The next section presents in more detail our main results and the organization of 
the paper. 
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2. Definitions and main results 



A planar map is a proper embedding of a connected graph into the sphere, where proper 
means that edges are smooth simple arcs which meet only at their endpoints. Two planar 
maps are identical if one of them can be mapped onto the other by a homeomorphism of the 
sphere preserving the orientation. A planar map is made of vertices, edges and faces. The 
degree of a vertex or face is the number of edges incident to it (counted with multiplicity). 
Following Bousquet-Melou and Schaeffer (2000), we consider a particular class of planar maps 
called constellations, see Fig[T] 



Definition 1. For p e [2, oo[, a (p-jconstellation is a planar map whose faces are colored black 
or white in such a way that : 

• adjacent faces have opposite colors, 

• the degree of any black face is p, 

• the degree of any white face is a multiple of p. 



Each edge of a constellation receives a canonical orientation by requiring that the white face 
is on its right. It is easily seen that the length of each oriented cycle is a multiple of p, and that 
any two vertices are accessible from one another. A constellation is rooted if one of its edges is 
distinguished. The white root face and black root face are respectively the white and black faces 
incident to the root edge. The root degree is the degree of the white root face. A constellation is 
said to be pointed if it has a distinguished vertex. Note that, in a pointed rooted constellation, 
the pointed vertex is not necessarily incident to the root edge. 

In this paper, we consider p-constellations subject to a control on white face degrees, i.e. 
for each positive integer k we fix the number of white faces of degree kp. This amounts to 
considering multivariate generating functions of constellations depending on an infinite sequence 
of variables (xk)k^i, where Xk is the weight per white face of degree kp and the global weight 
of a constellation is the product of the weights of its white faces. Two families of constellation 
generating functions will be of interest here. 

The first family is that of rooted constellations with a prescribed root degree. More precisely, 
for n > 1, let F n = F n (xi,X2, ■ ■ •) denote the generating function of rooted p-constellations with 
root degree np. By convention, we do not attach a weight x n to the white foot face and we set 
F = 1. 

The second family is that of pointed rooted constellations with a bounded "distance" between 
the root edge and the pointed vertex. More precisely, in a pointed constellation, we say that a 
vertex v is of type j e N if j is the minimal length of an oriented path from v to the pointed 
vertex (due to the orientation, it is slightly inappropriate to think of j as distance) and we say 
that an edge is of type j — > j' if its origin and endpoint are respectively of type j and j'. Then, 
for i 5= 1, let Vi = Vi(xi, X2, ■ ■ ■) denote the generating function of pointed rooted p-constellations 
where the root edge is of type j ' — > j — 1 with j ^ i. Note that V± is the generating function for 
rooted constellations (since for i = 1, the pointed vertex must be the endpoint of the root edge). 
Furthermore, we add a conventional term 1 to Vi, for all i > 1. 

The fundamental observation of this paper may be stated as: 
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Figure 1. A rooted 4-constellation of F 2 , with weight x\x% and the 4-excursion 
associated to it. 



Theorem 1. The sequences (-F n )n>o and {Vi)i^\ are related via the multicontinued fraction 
expansion 



(i) £ p,r 



1 



P- 1 \r. 
nStO 1 . T~~r v ii 



i-*n 



p-1 



i-*n 



ii=l T r Mi+fa 



P- 1 y. . . 

i 1 : i ^ j~~r y »i+»2+»3 



In the case p = 2, the r.h.s. of ([I]) reduces to an ordinary continued fraction (of Stieljes- 
type). This corresponds to the bipartite case discussed in (Bouttier and Guitter |2011 Eq. 



(1.13)): indeed, 2-constellations may be identified with bipartite planar maps upon "collapsing" 
the bivalent black faces into non-oriented edges. 

In the spirit of the combinatorial theory of continued fractions initiated by Flajolet ( |1980 1, 



an alternate formulation of Theorem [T] may be given in terms of some lattice paths. We call 
p-path a lattice path on Z x N made of two types of steps: rises (l,p — 1) and falls (1, — 1). 
Note that a p-path starting from (io, jo) only visits vertices (i, j) with i + j = i + j mod p. A 
p-excursion of length np is a p-path that starts at (0, 0) and ends at (np, 0) . It is well-known that 
such p-excursions are in one-to-one correspondence with p-ary rooted plane trees with n nodes, 
in number 1 +1 (" P T ^ 1 ) ■ To each fall in a p-path, we attach a weight Vi where i is the starting 
height of the fall (i.e. the fall starts from ( j, i) for some j, and thus ends at (j + 1)). We 

define the weight of a p-path as the product of the weights of its falls. As discussed in Section 
[3] the alternate formulation of Theorem [l] is then: 

Theorem 2. For all n ^ 0, F n is equal to the sum of weights of all p-excursions of length np. 

We prove this theorem in Section|4] As an illustration, for p = 3 and n = 1, 2, 3, the equality 
reads 

Fx = ^iF 2 

F 2 = V 1 V 2 {V l V 2 + V2V3 + V 3 V 4 ) 
^ F 3 = V^iVfVi + 2V 1 ViV 3 + VfVf + 2V 1 V 2 V 3 V 4 + 2V 2 V 3 2 V 4 + 

v 3 2 vi + v 2 v 3 v A v 5 + v 3 viv 5 + v 3 v 4 v s v 6 ) 
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We are now interested in inverting the relation ([!]), i.e. expressing Vi in terms of the F„'s. 
For p = 2, this may be done using Hankel determinants (see below). However, as soon as p 5= 3, 
it is not difficult to see that knowing the sequence (F„)„;>o alone is not sufficient (for instance 
in ([2| it appears that we have "twice as many unknowns as equations"). 

We thus need some extra knowledge. A possible way is to consider, for all n > and r 5= 0, 
the sum of weights of all p-paths that start at (— r, r) and end at (np, 0), which we denote by 

( r \ (Q\ ]\ 

Fn ■ Note that FA — F n and, since a p-excursion necessarily starts with a rise, F„ = F n +i- 
From now on, we restrict the values of r to the interval [0,p — 2], As will be discussed in Section 
those Fn ^ have a natural interpretation in terms of constellations. Let us write down the first 
3: 



few F^ for p 



(3) 



(i) 



VxV 2 {Vx + V 3 ) 



F 2 (1) = ViViQftVi + 2V 1 V 2 V 3 + V 2 V£ + VxVzVa + V^V± + V 3 V A V 5 ) 

It appears that, interlacing the equations in Q and ([2]), we obtain a triangular system of 
equations for Vi, .., V§. This fact turns out to be general and, furthermore, the solution to our 
inverse problem is provided by the following formula, which we derive combinatorially in Section 

El 

Theorem 3. For me [0,p — 1] and n e N ? we have the determinantal identity 

n ip-\-m 



(4) 



det F^ + i 



nru- 



where, for k e N, qk and r^ denote respectively the quotient and the remainder in the Euclidean 
division of k by p — 1, namely qk = Lj^rf J an d r k = k — (p — 1)<3%. 

Corollary 4. For me [0,p — 1] and n e N, we have 



(5) 



V m 



Hm\n— 1 ^m- l;r 
■ffo;n-Hp-l;n-2 



if m ^ 1 
if m = and n ^ 1. 



with the convention H m -\ = 1. 



Note that, for p 
m = or 1. For p = 



(6) 



= 2, we recover the Hankel determinants: 
3, the first few determinants read 



Ho-,0 = 




1 




Hi-fi = 


p(i) 


Hq-i = 




F (o) 

F (l) 




HiH = 


F (Q) 


Hq-2 = 


F (o) 


F (0) 


*f 
F« 
Ff 







,(0) 



det. 



i+j+m 



with 



2:0 



F. 



(0) 



#2:1 



F 



(0) 

i 



F v ; F 



r (0) 
2 

(1) 



and the reader may check, using ([2]) and ([3]), that those determinants are indeed monomials in 
the Vi's, in agreement with Q. We conjectu re that there i s a sim ple relation between H m - n and 



u< fn+i) p + m as defined at Equation (6.3f) in 
amounts to Corollary |4j 



Di Francesco 



(20051, so that Equation (6.30) ibid. 
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In this paper, we present a simple application of Theorem [3] to the case of Eulerian triangula- 
tions, namely 3-constellations where all the white faces are also triangular. This corresponds to 
specializing the generating functions defined above at the values p = 3, x\ = x, Xk = for k ^ 2, 
where a; is a variable controlling the number of triangles. In this case, the sequence (Vi)i^i is 
known to satisfy the simple recurrence 

(7) Vi=l + xVi (V 4 _i + V l+1 ) , i > 1, V = 0, 



sec 1 



Bouttier et al. (2003b). Clearly, this equation fully determines Vi as a power series in x. In 



the same reference, it was observed that the solution of this equation has a remarkably simple 
form 

(R) V = V (W)(W +4 ) 

[ ' 1 (l-^+l)(l-^+3)- 

where V and y are power series determined by the equations V = l + 2xV 2 , y+y -1 = {xV)~ l — 2. 
So far, there was no combinatorial explanation for the form of Equation [8] Theorem [3] provides 
such an explanation. Let us introduce the Fibonacci polynomials defined as follows (Flajolet 
|and Sedgewick] [20091 ec l- ( 62 )^ P- 327 ) 

(9) (p n+2 (z) = ip n+1 (z) - zip n (z), (p (z) = 0, (pi(z) = l, 

which are reciprocals of Chebyshev polynomials of the second kind, i.e. they satisfy the relation 

1 \ 1 - y n 



(10) LPn \y + y^ + 2) (l-y)(l + y)n-l- 
Then, Equation [8] results from the following proposition, which we prove in Section [5] 

Proposition 5. Let p = 3, x^ = x8^i for k 5= 1 and V = Xm=o n\(n+i)\ ( 2x ) w - F° r a ^ n ^ 0, 
we have 

Hojn-l = V^ 3n -^ 2 ip 3n+1 (xV) 

(11) Hi jn -i = V n ^ n+1 ^ 2 ^ n+2 {xV) 

(1 - xV)H 2 , n ^ = V^ 3n+ ^ 2 ^ 3n+3 (xV). 



3. p-PATHS AND MULTICONTINUED FRACTIONS 

We consider multivariate generating functions for p-paths with, for all i > 1, a weight Vi per 
fall starting from a height i. All results stated in this section will hold with (Vi)i^i a sequence of 
formal variables (i.e. its definition of Section [2] in terms of constellations is not needed). Recall 
that, for n, r > 0, F„ is defined as the generating function of p-paths starting from (— r, r) and 
ending at (np, 0). Let us add an extra weight t per rise and sum over all lengths, to obtain the 
generating functions 

(12) F^(t;V u V 2l ...):=J]F^t n . 



By elementary recursive decompositions of p-paths, we obtain the recursive equations 

'l + tF(P- 1 \t;V 1 ,V 2 ,...) ifr = Q 

V r F<°>(t; V r+1 ,V r+2 , . . .) FWfa V lt V 2 , . . .) if r > 1. 



(13) F^(t;V u V 2 ,...) 
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(For r = we remove the first rise. For r > 1 we perform first-passage decomposition at height 
j — 1.) We easily deduce that 

(14) F (0) (t; V 1 ,V 2 ,...) = — J- 

1 - tllLl Vi F n(t; v i+ i, v i+2 , . . .) 

By iterating this relation, we find that F^°\t; V\, V 2 , ■ ■ ■) is equal to the multicontinued fraction 
on the r.h.s. of ([!]). This shows that Theorem [2] (established in Section [4} implies Theorem [I] 




Figure 2. The unique configurations of non-intersecting lattice paths con- 
tributing respectively to H . 3 , £f 1;3 and H 2 -z for p = 3. 



We now turn to the proof of Theorem [31 which is an application of the celebrated LGV 



iotj 



lemma (Lindstrom 119691 Gessel and Viennotl 1989). We consider the weighted acyclic directed 



planar graph whose vertices are the (fc, I) e Z x N such that k + I = mod p, and whose edges 
are the rises (k,l) — > (k + 1, 1 + p — 1), weighted 1, and the falls (k, I) — * (k + 1, 1 — 1), weighted 
V;. For n e N, let and r n be respectively the quotient and the remainder in the Euclidean 
division of k by p — 1, and let 

(15) ^4„ := (-n - <?„,r„) = - r n ,r n ), B n := (np,0). 

The generating function for paths from A; to £?j (z, j 5= 0) is nothing but ^ij+V 

Theorem |3] then results from the LGV lemma and the following proposition, illustrated on 
Figi 

Proposition 6. Form e [0,p— 1], i/iere is a unique configuration of non-intersecting lattice 
paths connecting the sources A m , A m+ i, . . . , A m+n to the sinks Bq, B±, . . . , B n : for i e [0, 71], 
the source A m+ i is connected to the source Bi via the highest possible path, passing through 
(—to, to + ip) . The weight of this configuration is 0™=o Ilj^i" 1 • 

The proof is left to the reader. For p = 2 we recover the known combinatorial interpretation 



of Hankel determinants, see e.g. Viennot (1998]) 



4. Constellations, p-paths, and the slice decomposition 
In this section, we establish Theorem [2] and related results. 

Let us start with some definitions and notations. For n > 1, let !F n (resp. J 7 *) be the set 
of rooted (resp. pointed rooted) p-constellations with root degree np. The generating function 
of T n is F n as defined in Section [2j Similarly, let Vj (i > 1) be the set of pointed rooted p- 
constellations whose root edge is of type j — * j — 1 with j < i, to which we add a conventional 
"empty map" with weight 1. The generating function of Vi is V$. 

A constellated path is a p-path such that, for all i > 1, an element of Vi is associated with 
each fall from height i. The weight of a constellated path is the product of the weights of its 
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m+2 




Vs + l 



weight = 1 



weight = Vj 



Figure 3. The decomposition of a rooted pointed 4-constellation along left- 
most minimal path into slices (a). Slice corresponding to edge (m, m + 3) are 
reduced to a black edge (b) . Identifying the two sides of the boundary of a slice 
of an edge (m, m — 1) yields to a pointed-rooted constellation (c) 



associated constellations. We shall consider constellated excursions and bridges, where ap-bridge 
of length np is a p-path starting from (0, j) and ending at (np,j) from some j e N. Theorem [2] 
then follows from: 

Proposition 7. There is a weight-preserving bijection between T n and the set of constellated 
excursions of length np. More generally, there is a weight-preserving bijection between J 7 ' x N 
and the set of constellated bridges of length np. 

This bijection is a byproduct of the correspondence between constellations and mobiles in- 
troduced in Bouttier et al. (2004 1. Here we provide a direct construction, the so-called slice 



decomposition, extending the one given for p = 2 by |Bouttier and Guitter (2011). 

We start from a constellation C e J 7 *, and a £ e N (whose role will be discussed at the very 
end). The white root face of C forms a directed cycle which we denote by (vq, v\ , . . . , v np -\, v np = 
vq) where, say, "Uo is the endpoint of the root edge (hence v np -i is its origin). 

Let j s be the type of v s (recall that the type of a vertex v is the minimal length of an oriented 
path from v to the pointed vertex), and let 

(16) P:= ((Q,3o + e),(l,3i+e),...,(np, Jnp + l)). 
We claim that P is a p-bridge i.e. for all s e [0, np — 1], 

(17) j s+1 -j s e{-l,p-l}. 

Indeed, since (v s ,v s+ i) is an oriented edge, we have j s ^ j s+ \ + 1 by the definition of type. 
Furthermore, the black face incident to (v s ,v s+ i) forms an oriented path of length p — 1 from 
v s+ i to v s , hence j s+ i < j s + p — 1. Finally, since p divides the length of any oriented cycle on 
/, we have j s = j s+1 + 1 mod p, which concludes the proof of our claim. In particular, when C 
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is a rooted constellation that we point at the endpoint of the root edge, and £ — 0, then P is a 
p-excursion. 

Let us now explain how P is constellated. Assume that the constellation C is embedded in 
the plane with the white root face as outer face, and consider, for each s, the leftmost minimal 
(oriented) path from v s to the pointed vertex, where minimal refers to the length of the path. 
This family of paths leads to a decomposition of C into connected components that we call slices, 
see Figj3|a). A slice is associated with each edge (v s , u s +i) incident to the white root face (thus 
with each step of P). More precisely, it is delimited by (v s ,v s +i) and the two paths starting 
from v s and v s +i (in general, those two paths merge before reaching the pointed vertex, and we 
remove their common part). 

Observe that when j s +i = j s + p — 1, the slice is reduced to a single black face (as the path 
starting from v s+ i passes through v s after circumventing this face) with weight 1, see Fig. [3](b). 

Let us now assume that j s +i = j s — 1: we claim that the slice corresponds to an element of 
Vj s+ e- Two cases may occur. If the path starting from v s passes through v s +i, then the slice 
is empty, and corresponds to the empty map in Vj g+ e. Otherwise, let w s be the vertex where 
the two paths starting from v s and v s+ i merge. The boundary of the slice is therefore made of 
two (non-empty) oriented paths from v s to w s that do not meet except at their extremities. By 
construction those two paths have the same length: we may identify pairwise their edges, sec 
Fig|3](c). This identification preserve the degrees of the white faces, hence the weights, and the 
orientations guarantee that the resulting map is a constellation C s , which we root at (v s ,v s +i) 
and point at w s . By construction, the identified paths form in C s the leftmost minimal path 
starting with the root edge and ending at the pointed vertex. Thus, in C s the root edge is of type 
j' s —* j' s — 1, with j' s j s j s + £, so that C s e Vj s +e as claimed. Note that since the pointed 
vertex of C is incident to at least one white face, there is at least one s such that j' s = j s i.e. 
C s e Vj B \Vj s -i- This remark allows to characterize £ as the largest integer such that P — (0,1) 
is still a constellated bridge. 

In conclusion, the slice decomposition yields a constellated bridge. By following the steps in 
reverse order, it is clear how to reconstruct a pointed rooted constellation from a constellated 
bridge. By the above remark we also recover the integer £ and we may check that the corre- 
spondence is one-to-one (note that constellated bridges differing by a height shift yield the same 
constellation, but different values of £). Thereby we prove Proposition [7] and Theorem [2] 

Let us now mention other byproducts of our construction. 

Proposition 8. For all n ^ and re [0,p — 1], is the generating function of rooted 

p- constellations with root degree (n + l)p such that, if we denote by (vq, V\, . . . , Vr n+ i\ p = Vq) the 
directed cycle corresponding to the white root face (vq being the endpoint of the root edge), then 
all the vertices V\, ■ ■ ■ , w p _i_ r are bivalent. 

Proof. Apply the slice decomposition with vq as pointed vertex and £ = 0. The leftmost minimal 
path from v\ to vq visits successively «2, • • • ,v p - r , so that the constellated bridge starts with 
one rise followed by p — 1 — r falls corresponding to empty slices. Removing those p — r trivial 
steps, we obtain a constellated path with np + r steps going from height r to height 0. □ 



Proposition 9 (Bouttier et al. (20041; Di Francesco (20051). Let Fn be the generating 



function for constellated paths from (0, i — 1) to (np — We have: 

00 

(18) ^ = 1+^}] x n F^ 

n=l 

Proof. Apply the slice decomposition to a (non-empty) map in Vi- Taking £ = i — j, where 
j —* j — 1 is the type of the root edge, we obtain a constellated bridge of arbitrary length, 
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starting at height i — 1 and ending with a fall from height i, thus the factor Vi- The weight x n 
accounts for the white root face. □ 



Note the recursive nature of ( 18 ): F n l is a polynomial in the Vj's such that < (p— l)n. 
In particular, for i — > oo, Vi converges (in a suitable topology) to V the generating function of 
V = Ui>i an d i*n* -1 '*' — » (" 3 ^~ 1 ) yO 5-1 )- 1 , Hence V satisfies the recursive equation 

(19) - 1 ■ ^ ^P" 1 



r 



fc=l 



^ x n V n ( p ~ l \ 



By the Lagrange inversion formula, we may obtain an explicit expression for the coefficients of 
V (Bouttier 2005 p. 131), and recover the number of rooted constellations with a prescribed 
degree distribution (Bousquet-Melou and Schaeffer 2000 Theorem 2.3). 

By extending the construction in (Bouttier and Guitter 2011 Section 3.3), we may express 
the generating function F n in terms of V as follows. 

Proposition 10. For n 5= 0, and V the solution of (19), we have 

1 



(20) F n 



np + 1 



np + 1 

n 



V 



n(p-l) + l 



_y (y + l ( n P + 1 >. ( k P ~ lS \\ 

& \£ np + 1 \n-j J \k + jj) Xk 



y(k+n)(p- 



It would be interesting to have a similar formula for all F n r \ We could then apply Theorem [3] 
and Corollary [4] and deduce a general expression for Vi. We may express F n r ^ in terms of F n 
using a Tutte-like decomposition (Bousquet-Melou and Jehanne 2006 Section 5.3), for instance 
for p = 3 and r = 1 we have 



(21) 



F« = J] Xl F n+l + £ FiFn- 



/>-! 



(r) 

Applying Proposition 10 yields in principle to an expression of Fn in terms of V . In practice 
it seems to quickly become intractable, except in the case of Eulerian triangulations on which 
we focus next. 

5. Application to Eulerian triangulations 

We now specialize our formulas to the case of Eulerian triangulations, i.e. p = 3 and x^ = 
xSk i. Note first that (18) reduces to ([7]), while (19) yields the simple equation V = 1 + 2xV 2 , 
hence V is as in Proposition [5] Furthermore, (20) and (21) become respectively: 



(22) 



F 



\=(p^ ■{l-2xV)-p { ?_ 1 -{xV))v 



2n+l 



F n ^ = xF n+1 + £ F,F n . 



i=0 



where we introduce the shorthand notation p n 



(r) 



r+1 (3n+r+l\ 



3n+r+l 



the number of 3-paths from 



(— r, r) to (0, 3n). By the relation p n ^ = pffl + p„]_i (whose path interpretation is obvious), we 
may rewrite 

(23) f„ ■ (//„'" -ci -.--r ]-//„" ■ un) r- 



'„= (pW'(l-xV)-pW.(xV)) 



Using this relation in the expression (22) for Fn \ and using some summation formulas for the 
Pn s (following from several classical path decompositions), we arrive at 

(24) ' /..Ml ,1 .1 . i ..I I T/2n+2 



^ = (^.(1-^-^.^)) V 2 
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We are now ready to substitute these formulas into the determinants H n - m (m = 0, 1, 2) so as 
to prove Proposition [5j We will once again use the LGV lemma, to give another non-intersecting 
lattice path interpretation (NILP) of the H n . m 's. Let us consider the same acyclic directed 
planar graph introduced in Section [3] for p = 3 (i.e. the graph on which the 3-paths live), but 



now with a uniform weight 1 on every edge (rise or fall). Let A n and B n be as in (151 (for 
p = 3), namely A n = (— [3n/2j, r n ) and B n = (3n, 0), and let us extend the graph by adding the 
vertices 

(25) A' n := (-(3n+l)/2,-l), n > 

and the special weighted edges 
(26) 

j e 2 „ := A' 2n —> A 2n weighted 1-xV \e 2n +i ■= A' 2n+1 —> A 2n+1 weighted (1-xV) 2 
|e 2 „ := A' 2n -> A 2n+1 weighted -xV {e 2n+ i := A' 2n+1 -> A 2n+2 weighted -xV 

We readily see that the resulting graph is planar and, up to a factor V 2 ^ + ^' +1 (resp. l/ 2 ( 4 +J+ 1 )) j 
that the generating function for paths from A' 2i (resp. A' 2i+1 ) to Bj (i,j 5= 0) is nothing but 

Fi + j (resp. F^j). We may get rid of those extra factors by defining 

/ 97 x rp _ H0\k-1 rp _ Hi-k-i _ J?2;fc-l(l - XV) 

(41) J-3k+l ■ yk(3k-l)/2 ! J 3/c + 2 ■ yk( 3 k+l)/2 ' J 3/c+3- yfe(3fc+3)/2 

and it follows from the LGV lemma that T 3k+1 (resp. T 3k+2 and T 3 fc +3 /(l — xV)) is the gen- 
erating function of NILP configurations from the sources A' , . . . , A' k _ 1 (resp. A' 1 ,...,A' k and 
A 2 , . . . , A' k+1 ) to the sinks Bo, • • • , -Bfc-i- By convention, empty configurations (for n = 0) re- 
ceive the weight 1 so that T\ = T 2 = 1 and T3 = 1 — xV. We claim that the T n 's satisfy the 
recurrence relation 

(28) T n+3 = (1 - xV)T n+1 - xVT n 

and this is sufficient to prove that T„ = ip n (xV), since the Fibonacci polynomials with z = xV 
satisfy the same recurrence and initial conditions, by We prove (28) via a case- by-case 



analysis, depending on the residue of n modulo 3. Here we only check the case n = 3fc + 1, 
and leave the others cases to the reader (beware the extra 1 — xV factor in T 3k+3 ). Consider 
a NILP configuration from A' , . . . ,A' k to B , . . . , Bk, as counted by T 3 k+4- We first easily see 
that, for all i, the path i (from A\ to B{) visits (0,3i), and all the steps afterward are falls. We 
now distinguish whether the step before (0, 3k) on the uppermost path is a rise or a fall. 

• If it is a fall, then all the preceding steps are rises, so that the path visits A k +i and 
passes through the special edge e^, thus has weight — xV. All the remaining paths 
form an unconstrained NILP configuration from the sources A' , . . . ,A' k _ 1 to the sinks 
Bq, . . . ,Bk—i, with generating function T 3k +i- The global contribution of this case is 



-xVT 3k+1 (see Fig 4(a) I 



If it is a rise, i.e. the uppermost path visits (— l,3fc — 1), then, by non-intersection, the 
path % must visit (—1, 3z — 1) for all i 5= 1, and the path is reduced to the special edge eo, 
thus has weight 1—xV. Remove the path 0, and for each path i 5= 1, replace its steps after 
the vertex (— 1, 3i — 1) by 3i — 1 falls (see Fig (5j): the resulting configuration is a NILP 
configuration from the sources A' 11 . . . , A' k to the sinks B , . . . , B k _i, with generating 
function T 3k+2 . Since the correspondence is one-to-one, the global contribution of this 
case is — xVT 3k+2 (see Fig |4(b)| . 



This shows that (28 1 holds with n = 3k + 1. Together with the cases n = 3k + 2 and n = 3k + 3, 



this concludes the proof of Proposition [5j 



CONSTELLATIONS AND MULTICONTINUED FRACTIONS 
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The step before 
(0,3fc) is a fall. 




Si B 2 B 3 

A 3 A 2 Ai A 

(a) If the last step before 
(0, 3fe) is a rise, the upper-most 
path is deleted. 



The step before 
(0, 3k) is a rise. 




e 3 , (1-xV) 




Bn B\ B 2 



/;„ B 1 B 2 B 3 

A 3 A 2 Ai A es.(-xV) a'a'a' 

(b) If the last step before (0, 3fc) is a fall, for each path, we replace 
the steps after the vertex (— l,3i — 1) by 3i — 1 falls 



Figure 4. The decomposition of a configuration of i?o:4 into a configuration of 



Hq-3 (see (a) I or of Hi ;3 (see (b) ) 
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